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Using a tight binding model we calculate the conductance of monovalent atomic chains for 
different contact geometries. The leads connected to the chains are modelled as semi-infinite 
fee lattices with different orientations and couplings. Our aim is twofold: To check the validity 
of a three-parametric conductance formula for differently oriented leads, and to investigate the 
geometry dependence of the conductance oscillations. We show that the character of these 
oscillations depends strongly on the geometry of the chain-lead coupling. 
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1 Introduction 

Since Ohnishi p] et al. and Yanson et al. [2] demonstrated the existence of atomic chains, there 
has been a rapidly growing interest in this field [3-17]. The conductance of these chains was found 
to oscillate as a function of length, as it had been predicted earlier in Refs. [3] 0] and studied 
experimentally by e.g. Smit et al. [5]. This issue has also been investigated theoretically for alkali 
metals, and it was concluded that chains with odd numbers of atoms have a conductance equal 
to 2e 2 /h, while the conductance is smaller than 2e 2 /h for the chains with even numbers of atoms 
[HI [H [HI El] , (hereafter we refer to this as e < o). The opposite of this, namely e > o was found in 
noble metals [TO] and also in some calculations of alkali metals [I] . The role of symmetries was 
also studied in Refs. [Ill I12j and the authors stated that in chains with mirror symmetry and an 
odd numbers of atoms the conductance is exactly 2e 2 /h. However, Havu et al. [T3] showed for 
sodium wires that the oscillations may depend on the geometry. Hirose et al. [14] also pointed 
out that the oscillations are crucially influenced by the chain-lead coupling. 

There exist different simple models that help to understand and interpret the numerical and 
experimental results concerning the conductance oscillations, such as in [6l [TUl HS1 HB] , but none 
of them is suitable to explain the transformation between e < o and e > o due to the geometry 
dependence. 

In our previous work [17] we verified with the ab-initio code SMEAGOL [19l [20j [21] that a 
three-parameter formula can describe surprisingly well the conductance of monovalent atomic 
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Fig. 1 The lead(A)-chain(C)-lead(B) system. The coupling atoms(D) between the chain and 
the leads are described by scattering matrices with reflection coefficients r' L and tr those appear 
in Eq. |T}. 



chains, although the electron-electron interaction could destroy its applicability. In units of 2e 2 /h 
the conductance as the function of length is given by 



where d is the length of the translational invariant segment of the chain, k is the wavenumber 
of the transmitting electron, <&lr = &l + $r and $l(Qr) and r' L (r^) are the phase shift and 
magnitudes of reflection amplitudes at the left (right) side of the chain. Eq. fl} can be derived 
in the framework of ballistic transport, using the Landauer-Biittiker formula [22], considering two 
scattering regions (namely the coupling atoms between the chain and the leads), one conducting 
channel in the wire and a translational invariant region within the atomic chain. The scheme of 
this simple model can be seen in Fig. |T]). In our previous paper we investigated only symmetrical 
leads with one chosen orientation on the two ends of the chain and we modified the arrangement of 
the coupling atoms. However, in the experiments a chosen lead orientation can not be guaranteed. 
Therefore it is relevant to study whether the model can be applied to different lead orientations 
and couplings. This is the aim of the present paper. 

In section [5] we describe the method we use to calculate the conductance of atomic chains 
placed between fee leads. We explain how to obtain the Green's function of the leads, the Green's 
function of the atomic chain and how to compute the conductance of the total system. The 
numerical results are presented in section [3] We verify the validity of Eq. (Q]) by fitting the results 
and obtaining the reflection coefficients and the phases for systems with various leads orientations. 
If the ballistic model can be applied, then the fit parameters are not independent. For example 
by using two different symmetrical systems aa and (3(3 (for the definitions see subsection 12. 3|) it 
is possible to predict the conductance of an asymmetrical system a(3 with Eq. ([1]). We compare 
this prediction with the result of the direct numerical calculation for the a(3 system. We finish 
with the conclusions. 
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2 The model 

To calculate the differential conductance of an atomic chain placed between two leads we use a 
tight-binding model developed by Todorov et. al. [23]. The Hamiltonian, which has identical 
on-site energies and identical hopping energies reads 




if Rn Rfn ) 

(R„|i?|R rn ) = ^ 7, if R„ and R m are first neighbours, (2) 
otherwise. 

where |R„) is the atomic state of an electron located on the nth atom with position vector R„. 
Hopping is allowed only between the first neighbors and the topology of the atom arrangements 
determines the physical system. In the calculations Eq = and 7 = —1 is taken as the zero 
point and the unit of our energy scale. Three differently oriented blocks are investigated, i. e. the 
surface planes of the leads are (111), (100) and (110), so we have six different lead pairs. 

The two leads are formed from semi-infinite stack of planes specified by their normal vector. 
The chain and the two leads are described by the Hamiltonians H s and H 1 ' 2 , respectively, and 
coupled by the potential V. The retarded Green's function related to the total Hamiltonian H of 

the coupled systems is given by G(E) = \im e ^ + [E — H + ie\ . To calculate the conductance 
of the system, g, we start from the general expression derived by Todorov et al. [25] : 

g = - 1 -TT[p l {E F )i\E F )p 2 {E F )i{E F )}, (3) 

where p\ : 2 = — — ImG 1,2 are the density of states for lead 1 and 2 expressed via the corresponding 
Green's function G 1,2 and the transfer operator i is defined by i(E) = V + VG(E)V . To calculate 
the conductance we need then to determine the Green's functions G 1 , G 2 and G for the leads and 
for the coupled system. The total Green's function G can be obtained from the Dyson equation 



G = Go + G VG, (4) 



where 



G = G 1 +G 2 + G S , (5) 

and G s is the Green's function of the isolated chain. In the next two subsections we present the 
method we used to calculate the Green's functions of the leads G 1 ' 2 and the chain G s . 

2.1 The Green's function of the leads 

An fee crystal lead can be built up by stacking of its atomic planes. The Bloch state of an electron 
in one of such planes is given by 

m> = 4^E ex pN- R «)i R ™>' ( 6 ) 

* n 

where I denotes the Ith layer of the crystal, q is the transverse wave vector on the Ith layer and 
|R^) is the atomic state of the electron located on the nth atom of the Ith layer (having position 
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vector Rjj). Here 1/y/N is a normalization factor. The leads are supposed to be semi- infinite; 
I = 0, 1,2,.. .00. Both {|R4)} and {\l, q)} are complete orthonormal bases. Every atom has 12 
nearest neighbors. In the case of (100) planes there are 4 neighbors on the same plane while (2 x 4) 
on the next two adjacent planes. Similarly, for (111) planes there are 6 neighbors located on the 
same plane and (2 x 3) on the adjacent planes. When the lead is built from (110) planes there 
are 2 neighbors on the same plane, and (2 x 4) on the nearest and (2 x 1) on the second nearest 
adjacent planes. One can calculate the dispersion relation E(q) of the electrons on the layer as 
well as the coupling between layers, i.e., the matrix elements (I, q\H\m, q), where m = I, Z± 1, i±2 
depending on the orientation of the planes in the leads. To calculate the density of states that 
appears in Eq. we need to construct the Green's functions G 12 of the leads from the solution of 
the Dyson equation for the semi-infinite stack of the coupled planes. In the case of the orientation 
(111), the Green's function has already been derived by Todorov et al. [23]. In a similar way one 
can derive an explicit expression for the Green's function in the case of stacking of (100) planes. 
However, for stacking of (110) planes the method used in Ref. [23] cannot be applied because 
every atom has neighbors not only on the first but also on the second adjacent layers. In this case 
we used Haydock recursion [M] to calculate the Green's function numerically. The smoothing of 
the density of states is achieved by the commonly used termination of the recursion. The chain is 
connected to the leads in such a way, that the position of the first atom of the chain is equal to 
the atomic position in the next imaginary atomic plain of the lead. (Index I of this fictional plain 
would be equal to -1.) In the (111), (100) and (110) orientations there are 3, 4 and 5 connecting 
atoms, respectively. 

2.2 The Green's function of the atomic chain 

We now consider a one-dimensional atomic chain with n atoms labelled by S\,S2, • • ■ ,s n . In the 
tight-binding approximation the Hamiltonian of the chain can be written as 

n n—X 

H S = E J2 l*i><*i|+7E (h)( Sl+ i\ + h + i)( Sl \), (7) 

i=l i=l 

where are the orthonormal atomic basis in the chain. Due to the special form of H s , the 

Green's function G s (E) = lim £ _ f0 + LE — H s + ie^j can easily be calculated using recursively 
the following identity: 

This identity is valid for arbitrary vectors \x) and \y). This method is more versatile than the 
one used in Ref. [53]. It can be used, e. g. for atomic chains with dangling segments or loops of 
atoms. 



2.3 The conductance 

The conductance given by Eq. ([3]) can be obtained by performing the trace in the atomic basis 
{|R n )}. After straightforward manipulation we find 

9a,0 = ^ Im (/„) Im (fp) I ( Sl \G(E F )\s n ) \ 2 | 7 | 4 , (9) 

where a and j3 denote the orientation of the leads 1 and 2, i. e. (100), (110) or (111) planes. 
Factors f a and fp in Eq. ([9]) are related to the non-zero matrix elements of the density of states 
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of the leads. For the three different orientations we find the following explicit forms: 
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where G" m are the matrix elements of the Green's function between layers I and m of the 

a oriented lead. Finally, in Eq. @ (si\G(Ep)\s n ) are the matrix elements of the total Green's 
function of the system between the first and the last atoms of the chain. Using the Dyson equation 
(HJ we obtain 



( Sl \G(E F )\s n ) 



Gf r 



1 - (G{J a - G' m fp)\^ + D(G")f a f p \^' 



(11) 



where D(G S ) = (Gf 1 G^„ — Gf„Gf„) with matrix elements Gfj = (si\G s \sj) of the Green's func- 
tion for the atomic chain. We now have all the necessary formulas to calculate numerically the 
conductance. 

3 Numerical results 

The dispersion relation of the chain in tight-binding approximation reads 

E(k) = E + 2jcos(ka), (12) 

where a is the lattice constant. Thus, the chain has states only in the Eq — 2"f . . . Eq + 2j interval 
therefore this is the energy range where the system can work as a conductor. We have taken 
9 different energy points in this interval and calculated the conductance with 6 different lead 
pairs. As a first step we checked the validity of the ballistic transport model by doing some tests. 
We fitted Eq. (JTJ) on the numerical results to check whether the required relations between the 
parameters are fulfilled or not. For the phase $lr in Eq. {l} there holds an addition rule. The 
phase of the conductance oscillations of a chain between a and (3 leads is related to the phases of 
the oscillations between the symmetrical (aa and (3(3) leads as 



(13) 



This relation is fulfilled by our numerical results with an accuracy better than 1%. The values 
of the phases from the fit of Eq. (Q} are summarized in Table [TJ The good agreement is clearly 
visible at every energy level. The energy dependence of the phases can be seen in Fig. [2l For 
symmetrical leads r' L = tr, so these parameters can be determined also from the fit as it can be 
seen in Table Q] and in Fig. [3] 

Another proof of the validity of Eq. ([TJ) can be obtained as follows: By calculating the con- 
ductance of two different but symmetrical systems, for example between aa and (3(3 leads, with 
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Fig. 2 The energy dependence of the phases &lr for symmetrical leads. Signs +, x and o 
denote the orientations (111), (100) and (110), respectively. 



the fitted parameters, $ QCt , $/3/3, r a , rp, the conductance of the asymmetric a(3 system can be pre- 
dicted. To do this we need to substitute $ Q( g from Eq. (fl~3]) . r' L — r a and = rp in Eq. (JTJ . The 
comparison of the predicted and the directly calculated results is plotted in Fig. [4] and implies 
that differently oriented leads do not modify the availability of Eq. (JTJ) . The orientation of the 
leads influences the conductance essentially due to the different chain-lead coupling and different 
density of the atoms on the surface. The above results suggest that all of these properties can be 
described by the parameters r' L , rp and <&lr, and we do not need to modify the simple ballistic 
transport model for different leads. 

We found that the wavelength of the oscillations does not depend on the length and the character 
of the leads. It is essentially determined by the dispersion relation of the chain, similarly to our 
earlier results 17J. At half filled band the wavenumber is tt / '(2a), which gives two-atom-periodicity, 
but at other Fermi levels it becomes different. For example, for wavenumbers k = 7r/(4a) or 
k = (37r)/(4a) it yields four-atom-periodic oscillations, in agreement with the results of Thygesen 



Table 1 The phase $ and reflection amplitude r of the oscillations at different energies, resulted 
from the fit of Eq. ©. 



m 


$111,111 


$100,100 


$110,110 


$111,100 


$111,110 


$100,110 


riii, in 


^100,100 


^110, 110 


-1.96 


-2.976 


-2.741 


-2.591 


-2.861 


-2.781 


-2.661 


0.651 


0.706 


0.712 


-1.5 


-2.386 


-1.551 


-1.047 


-1.968 


-1.716 


-1.299 


0.199 


0.314 


0.368 


-1.02 


-1.300 


-0.626 


-0.151 


-0.972 


-0.745 


-0.391 


0.085 


0.233 


0.311 


-0.92 


-0.955 


-0.476 


0.024 


-0.696 


-0.446 


-0.226 


0.071 


0.224 


0.309 


-0.5 


1.101 


0.227 


0.547 


0.660 


0.810 


0.395 


0.067 


0.206 


0.307 


0.02 


2.336 


1.004 


1.093 


1.667 


1.718 


1.042 


0.114 


0.213 


0.323 


0.5 


2.787 


1.605 


1.515 


2.200 


2.140 


1.560 


0.177 


0.244 


0.353 


1.02 


3.033 


2.090 


1.985 


2.615 


2.440 


2.070 


0.272 


0.315 


0.415 


1.6 


3.124 


2.629 


2.434 


2.869 


2.779 


2.544 


0.456 


0.479 


0.559 
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Fig. 3 The energy dependence of the reflection coefficients r for symmetrical leads. Signs +, x and o 
denote the orientations (111), (100) and (110), respectively. 



8 « ® ® . 
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Fig. 4 The conductance (in units of 2e 2 /h) of a chain between asymmetric (a/3) leads, at 
energy Ef = O.O27 determined by Eq. ([T}, substituting the parameters of the symmetrical aa 
and f3f3 systems (crosses), and that from direct calculation (open circles) as a function of the 
length d (in units of the lattice constant a) for lead orientations a — (110) and ft = (100). 




5 10 15 20 25 5 10 15 20 25 5 10 15 20 25 

d d d 

a b c 

Fig. 5 The conductance (in units of 2e 2 //i) as a function of the length d (in units of the lattice constant 
a) of the chain between symmetrical leads at Ef — O.O27. Figs, a, b and c show the results for the lead 
orientations 111-111, 100-100 and 110-110, respectively. 



and Jacobsen |25j , who found that the three conductance channels of aluminium have wavenumbers 
7r/(4o), 7r/(2o), and (37r)/(4a), and the conductance oscillations have four atoms periodicity. The 
periodicity of the oscillations given by Eq. ([1]) is also in agreement with the earlier prediction of de 
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la Vega and co-workers [55] , who predicted that every conduction channel oscillates as a function 
of length with a wavenumber 2k p. 

For monovalent atomic chains the charge neutrality condition requires wavenumber kp = 
tt /(2a). Then, it follows from Eq. (p} that the conductance is an alternating function of the 
number of atoms. At hp — n/(2a) the phase can not be determined easily from the fit on nu- 
merical data and it seems to be an unnecessary parameter. Maybe this fact can explain why the 
phase was missed in earlier interpretation models [5J [TUl HH HE]. However, in spite of these 
difficulties the phase is very useful to interpret the geometry dependence of the oscillations. To 
avoid numerical problems, we could take the limits lmiE^E F and lim^^^p r(E), but as 

can be seen in Figs. [2] and [3] these parameters are smooth functions of the energy and therefore 
we choose Ep — O.O27 in the calculations. For symmetrical leads some numerical results are 
plotted in Fig. [5] Because of the different groups of atoms connecting the chain to the lead, the 
oscillations have different phases. Between leads with planes (100) and (110) one can see that 
the conductance for chains with an even atoms is smaller than for chains with an odd number of 
atoms, i. e. e < o. However, in the case of (111) orientation, the so called even-odd effect changes 
its sign and e > o. This demonstrates the sensitive geometry dependence of the conductance. 
The analysis of Eq. (Q]) implies that the maximum value of the conductance can be (2e 2 /h) only 
if the cosine term is equal to 1 and therefore the phase is 0. It is also clear from the figure that 
not only the phases, but also the amplitudes depend on the coupling, since according to Eq. |T]) 
larger reflection coefficient r yields larger amplitude of the oscillation. Indeed we have found the 
relation nu < rioo < ^no, (0.114, 0.213 and 0.323, respectively at Ep = 02027) which follows the 
same pattern as the amplitudes of the oscillations. In the case of orientations (111), (100), (110) 
there are 3, 4 and 5 coupling atoms between the lead and the chain, respectively, suggesting that 
more atoms may increase the reflection amplitudes. However, it is also necessary to stress that 
the different chain-lead couplings result in different phases, and these can modify the amplitudes 
of the oscillations also. 

4 Conclusions 

We have shown that the differently oriented leads do not affect the applicability of Eq. ([1]) . How- 
ever, the geometrical arrangement of the atoms connecting the chain to the lead appreciably 
influences the conductance, since the phases and the amplitudes of the oscillations change signifi- 
cantly and the so-called even-odd effect can change its sign. This implies that numerical results of 
different calculation methods are comparable only in the case of equal geometrical circumstances. 
The validity of Eq. can advance further investigations of the geometry dependence of the os- 
cillations, since it can be used to predict the transport characteristics of asymmetrical leads, so 
that it would only be necessary to perform simulations of symmetrical systems. 

Since the parameters of the oscillations depend so sensitively on the details of the coupling 
between the chain and the leads, it is not clear at first sight why the even-odd effect can be de- 
tected experimentally. A possible explanation is the following. In the experiments the average of 
conductances of several different contacts is measured. The different geometries and consequently 
oscillations with different phases could suppress the even-odd effect. However, there are energeti- 
cally preferred surface orientations for a chosen material, and these facets assist the formation of 
preferred chain-lead couplings, resulting favored certain r and $ parameter values. Rodrigues et 
al. [27] [28] investigated the effect of these preferred surfaces on the conductance of point contacts 
experimentally, and they found that in the case of gold the (111) facets are more favorable. An 
other argument can be the following. In the experiments after one breaking of the contact the tips 
are brought together again to form a new contact. The displacement and variation of both parts 
are usually small. Therefore the recreated contact has similar geometry and presumably similar 
oscillation parameters. 
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In summary we showed that Eq. ([T]) and specially its phase plays an important role to interpret 
the conductance oscillations of atomic wires of monovalent elements. 

Acknowledgements The authors would like to thank the helpful discussions with Gergely Szirmai and 
Tamas Herpay. VMGS thanks the European Union for a Marie Curie grant. 

References 

[I] Hideaki Ohnishi, Yukihito Kondo, Kunio Takayanagi, Nature 395 780 (1998) 

[2] A. I. Yanson, G. R. Bollinger, H. E. van den Brom, N. Agrai't, J. M. van Ruitenbeek, Nature 395, 
783 (1998); 

[3] P. L. Pernas, A. Martm-Rodero, and F. Flores, Phys. Rev. B 41, R8553 (1990). 
[4] N. D. Lang, Phys. Rev. Lett. 79, 1357 (1997). 

[5] R. H. M. Smit, C. Untiedt, G. Rubio-Bollinger, R. C. Segers, J. M. van Ruitenbeek, Phys. Rev. Lett. 
91, 076805 (2003). 

[6] H.-S. Sim, H.-W. Lee, and K. J. Chang, Phys. Rev. Lett. 87, 096803 (2001). 
[7] R. Gutierrez, F. Grossmann, and R. Schmidt, Acta Phys. Pol. B 32, 443 (2001). 
[8] S. Tsukamoto and K. Hirose, Phys. Rev. B 66, 161402(R) (2002). 
[9] Y. Egami, T. Ono, and K. Hirose, Phys. Rev. B 72, 125318 (2005). 

[10] Y. J. Lee, M. Brandbyge, M. J. Puska, J. Taylor, K. Stokbro, and R. M. Nieminen, Phys. Rev. B 69, 
125409 (2004). 

[II] H.-W. Lee and C.-S. Kim, Phys. Rev. B 63, 075306 (2001). 
[12] Z. Y. Zeng and F. Claro, Phys. Rev. B 65, 193405 (2002). 

[13] P. Havu, T. Torsti, M. J. Puska, R. M. Nieminen, Phys. Rev. B 66, 075401 (2002) 
[14] K. Hirose, N. Kobayashi, and M. Tsukada, Phys. Rev. B 69, 245412 (2004). 
[15] E. G. Emberly and G. Kirczenow, Phys. Rev. B 60, 6028 (1999). 
[16] H.-W. Lee, H.-S. Sim, D.-H. Kim, and K. J. Chang, Phys. Rev. B 68, 075424 (2003). 
[17] P. Major, V. M. Garcfa-Suarez, S. Sirichantaropass, J. Cserti, C. J. Lambert, J. Ferrer, G. Tichy 
Phys. Rev. B 73 045421 

[18] Recent experiments have also studied the interaction of atomic chains with simple molecules; Sz. 

Csonka, A. Halbritter and G. Mihaly, Phys. Rev. B 73, 075405 (2006); W. H. A. Thijssen, D. 

Marjenburgh, R. H. Bremmer, and J. M. van Ruitenbeek Phys. Rev. Lett. 96, 026806 (2006). 
[19] A. R. Rocha, V. M. Garcfa-Suarez, S. W. Bailey, C. J. Lambert, J. Ferrer, and S. Sanvito, Phys. 

Rev. B 73, 085414 (2006). 

[20] A. R. Rocha, V. M. Garcfa-Suarez, S. W. Bailey, C. J. Lambert, J. Ferrer and S. Sanvito, Nature 
Materials 4, 335 (2005). 

[21] V. M. Garfa-Suarez, A. R. Rocha, S. W. Bailey, C. J. Lambert, S. Sanvito, and J. Ferrer, Phys. Rev. 
Lett. 95, 256804 (2005). 

[22] R. Landauer, IBM J. Res. Dev. 1, 223 (1957); M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986); S. 

Datta, Electronic Transport in Mesoscopic Systems, (Cambridge University Press, Cambridge, 1995). 
[23] T. N. Todorov, G. A. D. Briggs, and A. P. Sutton J. Phys. Condens. Matter 5, 2389 (1993); T. N. 

Todorov, Phys. Rev. B 54, 5801 (1996). 
[24] R. Haydock, The recursive Solution of the Schrodinger equation in Solid State Physics edited by 

H. Ehrenreich, F.Seitz, and D. Turnbull, vol 35, pages 215-294 Academic press, New York, 1980; 

R. Haydock, V. Heine, and M. J. Kelly, J. Phys. C 5, 2845 (1972); [INSPEC] ibid. 8, 2591 (1975); 

[INSPEC]; C. Lanczos, J. Res. Natl. Bur. Stand. 45, 225 (1950). 
[25] K. S. Thygesen and K. W. Jacobsen, Phys. Rev. Lett. 91, 146801 (2003). 
[26] L. de la Vega, A. Martin Rodero, A. Levy Yeyati, A. Saul Phys. Rev. B 70, 113107 (2004) 
[27] V. Rodrigues, J. Bettini, A. R. Rocha, L. G. C. Rego, D. Ugarte, Phys. Rev. B 65, 153402 (2002) 
[28] V. Rodrigues, T. Fuhrer, D. Ugarte, Phys. Rev. Lett. 85 4124 (2000) 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Wcinhcim 



